Endoscopic Tomography and Quantum-Non-Demolition 

Mauro FortunatoQ and Paolo Tombesi 
Dipartimento di Matematica e Fisica, Universita di Camerino, 1-62032 Camerino, Italy 
and Istituto Nazionale per la Fisica della Materia, Unita di Ricerca di Camerino 



oo 
o 

Q 



> 

m 
o 

00 ■ 

o\ ■ 

^ : 

Ok 
^— > ' 
G ■ 
cd 

3 : 
cr 



Wolfgang P. Schleich 
Abteilung fur Quantenphysik, Universitat Ulm, D-89069 Ulm, Germany 
(Received February 1, 2008) 

We propose to measure the quantum state of a single mode of the radiation field in a cavity — 
the signal field — by coupling it via a quantum-non-demolition Hamiltonian to a meter field in a 
highly squeezed state. We show that quantum state tomography on the meter field using balanced 
homodyne detection provides full information about the signal state. We discuss the influence of 
measurement of the meter on the signal field. 

PACS numbers: 03.65.Bz, 42.50.Dv 



I. INTRODUCTION 

How to measure the quantum state of a single mode of 
the radiation field in a cavity? Various possibilities 0-0] 
offer themselves. However, a straightforward application 
of the method of quantum state tomography suggested in 
Rcf. J8| and implemented experimentally in Refs. |] Q 
does not work, since by coupling the field out of the res- 
onator we change the field state, fn the present paper we 
propose to couple the field via a quantum-non-demolition 
(QND) interaction |l2| to a meter field on which we then 
perform tomography using a balanced homodyne detec- 
tor. In this way we combine the idea of probing, that 
is doing endoscopy on the field without taking it out of 
the cavity, and the tool of tomography and arrive at the 
method of endoscopic quantum state tomography. 

The goal of the present paper is to obtain informa- 
tion about the full quantum state of a single mode of 
the radiation field. To bring out the physics most clearly 
we assume that this field, referred to in the remainder of 
this article by the signal mode, is in a pure quantum state 
and neglect damping. We emphasize, however, that the 
method presented here also applies to a signal field de- 
scribed by a density operator. In contrast to the method 
of quantum state tomography |§-|ll[ based on homodyne 
detection, the present technique does not couple the sig- 
nal field out of the resonator. In order to measure the 
signal field we couple it in a linear way to a meter field. 
Moreover, we couple both to a pump field. This allows 
us to achieve a quantum-non-demolition Hamiltonian de- 
scribing the interaction between the signal and the meter 
mode. The use of a QND-Hamiltonian suggests that one 
might be able to arrange the scheme in such a way as to 
measure a complete quadrature distribution without re- 
preparing the quantum state. In other words, repeated 
measurements on the meter change the signal state but 
keep the quadrature distribution invariant. We show that 



unfortunately this is not the case. This is closely related 
to the question if the wave function of a single quantum 
system could be measured Indeed Ref. [Q sug- 

gests that the wave function of a single quantum system 
could be measured by employing a series of "protective 
measurements" where an a priori knowledge of the wave 
function enables one to measure this wave function and 
protect it from changing at the same time. However, Al- 
ter and Yamamoto (l5) showed that a series of repeated 
weak quantum non-demolition measurements gives no in- 
formation about the wave function of the system. The 
same authors have also argued that it is not allowed 
to measure the full state of a single quantum system. 
Recently, D'Ariano and Yuen (T?]] have independently 
proven the impossibility of measuring the wave function 
of a single quantum system. The present intentions are 
much less ambitious since, eventually, we do not want to 
measure the full state of a single quantum system, but 
only the quadrature probability distribution. 

The article is organized as follows: In Sec. [n] we re- 
derive the relevant QND Hamiltonian emphasizing its 
dependence on the phase of the pump field which allows 
us to probe all quadratures of the signal field. We devote 
Sec. Ill to the calculation of the entangled state of meter 
and signal originating from the unitarytime evolution 
due to the QND Hamiltonian. In Sec. [TV] we study the 
influence of the measurement of the meter on the sig- 
nal field and in Sec. ^ we consider two special cases: in 



VI 



we then 



phase and out of phase measurements. In Sec. 
turn to the question of tomography using a QND Hamil- 
tonian. In Sec. VII we give a general argument which 
shows the impossibility of having a (QND) measurement 
which simultaneously keeps the probability distribution 
unchanged and gives information about the measured ob- 
servable. We conclude in Sec. VIII by summarizing our 
main results. In order to keep the article self-contained 
we have included all relevant calculations but have sum- 
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marized longer ones in Appendices ^ and [b]. 



II. QND HAMILTONIAN 

In the present section we derive the QND Hamiltonian 
used in our tomographic scheme to couple the signal to 
the meter field. This treatment brings out clearly how 
the phase of the pump field allows us to probe every 
quadrature of the signal. 

Our model starts from the Hamiltonian 

H = ihxiala^&p - a s a m a^\ + ifia[a\a m - a s a\„\ , (2.1) 

where a s {d\), a m (al n ), and a p {a/ p ) denote the annihila- 
tion (creation) operators of the signal, meter, and pump 
field, respectively. The parameters \ an d & measure the 
coupling between the three fields, and the meter and sig- 
nal field, respectively. 

A possible scheme of the measurement strategy sug- 
gested in this paper is shown in Fig. [lj We assume that 
the crystal is present in the cavity when we prepare the 
signal field. In this case the pump and the meter field 
are in vacuum states and the resulting modifications on 
the signal due to the presence of the crystal can be easily 
taken into account. 

When the pump field is highly excited we can describe 
it by a coherent state of amplitude a and phase 2(f>, that 
is 



(2.2) 



Here we have defined the phase 2<fi rather than 4> as to 
simplify the resulting equations. It is the variation of 
this phase <j) of the pump field which allows us to per- 
form tomography on the signal field. To understand this 
in more detail we substitute the coherent state approxi- 
mation, Eq. ( [2.2[ ), of the pump field into the Hamiltonian, 
Eq. (2.1), and find after minor algebra 



H = 2hax s {4> + tt/2) • x m {(j>) 



(2.3) 



Here we have arranged the strength a of the pump field 
such that x a = a - Moreover, we have introduced the 
quadrature operators 



Xj(0) 
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ate- 



(2.4) 



of the signal (J = s) and the meter (j = m) mode at the 
angle 9. 

Note that due to the special choice x a = c of the 
pump field we have achieved an interaction between the 
signal and the meter which couples the quadrature op- 
erator £ m (</>) of the meter at phase angle <p to the out- 
of-phase quadrature operator x s (4> + tt/2) of the signal. 
Such Hamiltonians have been studied extensively 18 24j 



in the context of quantum non-demolition measurements. 
In the present paper we analyze how such a Hamiltonian 
can be used to measure the quantum state of the signal 
field. We note that according to the QND Hamiltonian 
Eq. fl2.3| ) a measurement of the meter at a fixed phase (j> 
of the pump field provides information about the signal 
in the out of phase quadrature. By varying the phase <f> 
of the pump field we can probe in this way all quadra- 
tures of the signal. We conclude this section by noting 
that we can achieve a measurement of the meter quadra- 
ture operator by a homodyne measurement of the meter 
mode. 



III. ENTANGLEMENT 

We now calculate the combined state |^) of signal and 
mete r ob tained from the QND interaction Hamiltonian, 
Eq. @. 

When we couple the signal and meter mode prepared 
initially in the states \ip s ) and \ip m ) we find the quantum 
state 



|¥(t)> =exp(-iHt/h)\i/j m )\A) 

= ex.p{-2icrtx s (<f> + n/2)x m ((f))}\'ip m )\^ s ) (3.1) 

for the combined system after the interaction time t. This 
time is determined by the decay time of the cavity. 

To evaluate the above expression we expand the ini- 
tial signal state in quadrature states \x s {(j) + tt / 2)) of the 
phase angle 9 S = <f> + tt/2, that is 



\1>. 



dx s ip s (x a ;<p + ?r/2)\x s (<f> + Tr/2)) . (3.2) 



We emphasize that this representation and, in particular, 
the wave function ip s (x s ;<f) + tt/2) = (x s (4> + x/2)\ip a ) 
depend crucially on the angle 6 a . 

We substitute the expression Eq. ( |3 . 2| ) for the signal 
state into Eq. (fO), use the eigenvalue equation 



x s {6)\x s {6))=x s \x a {9)) 



(3.3) 



for the signal quadrature state \x s {9)) at angle 9, and 
arrive at the combined state 



!*(*)) 



dx s tp s (x s ; (j> + -k/2)\x s (4> + tt/2)) 



x exp[-2icrtx s x m ((j))]\ip r , 



(3.4) 



of signal and meter. 

To find the action of the exponential operator in 
Eq. (3.4) on the meter state \ip m ) h is convenient to ex- 
pand \ipjn) in quadrature states \x m (9)) of the meter at 
the angle 9, that is 
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(3.5) 



where ip m (x m ;6) = (x m (9)\ip m ) denotes the wave func- 
tion of the meter state at the angle 9. Note that this 
angle is still arbitrary and is not necessarily identical to 
the angle <fi in the Hamiltonian. According to the Ap- 
pendices |X] and ^ we find 



exp [-i(2(xtx s )x m (<P)]\i>v 



dx m exp[-i7(x s , x m ;9 - <f))] 
x tJj 2atx s sm(0 -0); 8}\x m (0)) 



(3.6) 



where 



(atx s ) 2 sin [2 (9 
+2crtx 

S&TYL COS! I 



(3.7) 



denotes the phase accumulated due to the interaction. 
Hence the combined quantum state reads 



l*(*)) 



dx. 



dx m ip s (x s ;4> + 7r/2) 



xipm[xm - 2atx s sin(6» -(f); 9] 

x exp[-vy(a; a , x m ; - (f>)} 

x\x s {4> + K/2))\x m {6)) . (3.8) 

We note that due to the coupling between the meter 



and the signal via the Hamiltonian Eq. (2.3), the meter 
wave function ip m (x m ; 9) at the angle 9 gets shifted by 
an amount Sx m = 2atx s sin(0 — 4>). This shift is propor- 
tional to the interaction strength erf, the signal variable 
x s and the sine of the angle 9 — <f>. 



IV. SIGNAL STATE CONDITIONED ON METER 
MEASUREMENT 

In the preceding sect ion we have calculated the entan- 
of the combined system. In the 



gled state |*), Eq. (|3 



present section we show how a measurement of the meter 
influences the state of the signal. In particular, we use 
the Wigner function approach to discuss the properties 
of the signal state conditioned on a quadrature measure- 
ment of the meter variable. Here we first consider an 
arbitrary quadrature state of phase angle 9 and then in 
Sec. |v| focus the discussion on two special cases. 
According to Eq. (|3.8|) the conditioned state 



l^ c) > 



f 



VW(x^j 



(x m (9)\V(t)) 



(4.1) 



of the signal given that our quadrature measurement at 
angle 9 has provided the value x m reads 

oo 

|^ c) ) = J dx s ip s (x s ;<t) + Tr/2)f(x s \x m )\x s (cl> + Tr/2)) , 

— OO 

(4.2) 



where the filter function 
1 



f{x s \ x Ti 



--ipm[x m ~ 2atx 8 sin( 



y/W(x m ) 
x esxp[-ij(x s , x m ; 9 - (/>)} 



(4.3) 



originates from the interaction of the signal with the me- 
ter. The probability W(x m ) of finding the meter variable 
x m follows from the normalization condition 



i = (v4 c) i^ c) >, 



(4.4) 



that is 



W(x m )= I dx s \ij s (x s ;<j> + Tr/2)\ 2 

x \ip m [x m ~ 2atx s sin(6 - 4>);9] | 2 . (4.5) 



Equation (4.2) clearly shows how the measurement of the 
meter influences the quantum state of the signal: The 
filter function determined by the wave function of the 
meter selects those parts of the signal wave function that 
are entangled with the corresponding parts in the meter. 
To study this in more detail we now calculate the Wigner 
function lEal 



oo 

W^ c) {x Sl p s \x m ) = i J dyt 



■y/2\1> 



(c)\ 



x(^\x s +y/2) (4.6) 

of the signal state conditioned on the measured meter 
value x m . For the sake of simplicity we have suppressed 
the angle (j> + ir/2 at the quadrature states \x + y/2) and 
\x- y/2). Substituting the state |v4 c) ), Eq. (p~^), into 
this expression we arrive at 

oo 

W s c \x s>Ps \x m ) = ±- I dyj™i> 8 {x s -y/2) 

— OO 

x ip*{x s + y/2)f(x s - y/2\x m )f*{x s + y/2\x m ) . (4.7) 
We express the integral as the convolution pq | 

W^(x s ,p s \x m ) = 

dp'W s (x s ,p s -p')W f (x s ,p'\x m ) (4.8) 
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between the Wigner function 

W s (x s ,p s ) = 

oc 

i- I dye^y^ s (x s ~y/2)r s ^s+y/2) (4.9) 



of the meter via the relation 
1 



Wf(x s ,p\x m ) 



W m [x m - 2atx s sin( 



W{x m ) 

p +2atx s cos(9 - (j>) + x m cot(6» - (j>)) . (4.15) 



of the original signal state and the Wigner function 

Wf(x s ,p s \x m ) = 

oo 

^ f dye^yf(x s ~y/2\x m )r(x s +y/2\x m ) (4.10) 



of the filter provide d by the meter measurement. We 
can e asily prov e Eq. ( |4.8| ) by s ubst ituting the expressions 
Eqs. ( |4.9| ) and ( 4.1C ) into Eq. (h8), interchanging the in- 
tegrations and performing one of them using the resulting 
delta function. We then indeed recover the integral (4.7). 

If we substitute the filter function Eq. (4.3) into the 
Wigner function Eq. (4.10), after minor algebra we obtain 



Wf(x s ,p s \x m ) 



1 



1 



W{x m ) 2tt 



dy 



x exp{iy[p s + 2(at) 2 x s sin[2(6> — </>)] + 2atx m cos(0 — 4>)]} 
xipm[x m - 2atsin(9 - 4>){x s - y/2)] 
xiP^[x m -2atsm(6-(j>)(x 8 +y/2)} . (4.11) 

When we introduce in the last integral the new integra- 
tion vari able y = 2at sm(9 — <fi)y, and in the convolution 
Eq. (4.8) the integration variable p = p' /[2at sin(9 — </))], 
the Wigner function of conditional state 

oo 

W ( s c) (x Sl p s \x m ) = I dpW s [x s ,p s - 2atsm(9 - 4>)p] 



xW f (x s ,p\x m ) 



(4.12) 



is the convolution of the Wigner function W s [Eq. (4.9)] 
of the signal state and the Wigner function 



Wf(x s , p \x m ) 



1 



W(x m ) 2tt 



dy 



x exp[iy(p + 2atx s cos(9 — <fi) + x m cot(9 — (/>))] 

x tp m (x m - 2atx s sin(6> - <fi) + y/2) 

x V40*m - 2atx s sm(9 - 4>) - y/2) (4.13) 

of the filter function. We express the latter in terms of 
the Wigner function 

W m (x m ,Pm) = 



2vr 



dye ip ™ y ^ m (x m - y/2)r m (x m + y/2) (4.14) 



V. SPECIAL EXAMPLES FOR CONDITIONED 
SIGNAL STATES 



Whereas in the discussion of Sec. IV the angle 9 of 



the meter quadrature is still arbritrary, we concentrate 
in the present section on two distinct cases: We choose (i) 
9 — <f>, that is we measure in phase and (ii) 9 = <p + tt/2, 
that is out of phase measurement. 



A. In phase measurement 



If we choose the angle 9 of the m eter quadrature to be 



identical to <fi, the state l^), Eq. (3 
system reduces to 



of the complete 



!*(*))= J dx s J dx m il] s (x s ;(f) + '7r/2)ip m (x m ;4>) 

— oo — oo 

xexp(— i2(xtx s x m )\x s {ct> + ^/2))\x m {<t>)) . (5.1) 

Here we have made use of the phase 7 = 2atx s x m , 
Eq. (3.7), for 9 = <\>. Note that this expression also fol- 



lows immediately from the Hamiltonian Eq. (2.3) and the 
expansions Eqs. (3.2) and (3.5) of the meter and signal 



states. We emphasize that in this case the meter wave 
function is not shifted. Nevertheless, th e tw o states are 
still entangled via the exponential, Eq. (3^). Since the 
shift &x m vanishes, the probability 



W(x m ) = |^ m | 2 / dx s \ip s {x a )\ 2 = IVJ 1 



(5.2) 



of finding the meter variable x m following from Eq. ( |4.5| ) 
for 9 — cf> is identical to the initial probability of the 
meter, that is 



W(x m ) = \ip m (x m )f 



(5.3) 



Here we have used the fact that the original signal 
wave function is normalized. Hence, up to an over- 
all phase fi m determined by the meter wave function 
^(xm) = ^(xm)! exp[i/j,(x m )], we find from Eq. (h3) the 
filter fun ction f(x s \x m ) = exp(— i2atx s x m ), and from 
Eq. (h2) the conditioned signal state 



00 



l^s ) — / dx s ips(x s ) exp(-i2atx s x m )\x s ((j> + 7r/2)) 



(5.4) 
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Note that the measurement of the meter has indeed 
changed the state of the system but did not alter the 
probability 

W(x s ) = |(x s |^ c) )| 2 = \^ s {x s )e^ V {-i2atx s x m )\ 2 

= \^s(x s )\ 2 (5.5) 

of finding the signal variable x s . This effect of the me- 
ter measurement comes out most clearly in the Wigner 



function W* of the conditioned system state, Eq. (4 



From Eq. (4.11) we realize that for 9 = <j) the Wigner 
function of the filter reduces to a delta function in the 
momentum shift, that is 

Wf(x s ,p s \x m ) = S(p s + 2atx m ) , (5.6) 

and the Wigner function following from the convolution 



Eq. (4.8) reads 



W^(x s ,p s \x m ) = W s (x s ,p s + 2atx m ) 



(5.7) 



Hence, the measurement has left untouched the shape of 
the original state represented here by the Wigner func- 
tion but has moved it along the momentum axis by an 
amount of 2atx m . Consequently, the measurement did 
not change the probability distribution in the conjugate 
variable, namely the x s variable. We note, however, that 
in this way we cannot gain information about the sig- 
nal since according to Eqs. ( |5.2D and (5.3) the probabil- 
ity distribution W(x rn ) of measuring the variable x m is 
identical to the original distribution. 

This finding is actually a rather general result. In fact, 
it can be rigorously shown [^7j that a (QND) measure- 
ment which does not change the probability density of the 
observable which is being measured on a single quantum 
system gives no information about the measured observ- 
able. Its proof, restricted for cl arity to the model consid- 
ered here, can be found in Sec. VII. 



B. Out of phase measurement 

We now turn to the case of 9 = <p+ir/2. In this case the 
shift Sx rn = 2atx s in the meter wave function is maximal 
and according to Eq. (|j]) the phase 7 vanishes. Hence, 
the combined state 

00 00 

|*(t))= J dx s J dx m iP s (x s ;<j) + TT/2) (5.8) 

— OO —OO 

2atx s ;cj> + tt/2) \x s {<j> + tt/2)) \x m ((j> + tt/2)) 

is an entangled state in which the entanglement between 
the meter and signal is due to the shift of the meter. 
In contrast to the discussion of Sec. VA we can now 



deduce properties of the signal from the shift of the meter 
wave function. Unfortunately, we cannot simultaneously 



keep the probability distribution W(x s ) — \i[) s (x s )\ 2 of 
the original signal state invariant, in accordance with the 
discussion at the end of Sec. VA (see Sec. |VII| ). Indeed, 
we find from Eqs. (4.2) or (5.4) the conditional state 



IV'. 



(c)\ _ 



1 



W(x n 



dx s ip s (x s )ijj m (x m - 2atx s ) \x s ) 



(5.9) 



of the system given the meter measurement at phase <fi + 
tt/2 has provided the value x m . The probability 



W(x m ) 



dx a \ip„(x a )\ 2 \ip m (x m - 2atx s )\ 2 (5.10) 



of finding the meter value x m following from Eq. (4.E) is 
now a convolution of the system and the meter function. 
In Sec. VI B we will use this relation to perform tomog- 
raphy on the system. However, in the present section we 
focus on how the measurement influences the signal state. 
We note that in contrast to the discussion of Sec. V A the 



meter measurement has changed the conditional distribu- 
tion 



W^(x s \x m ) = \(x s \^)\ 2 

\ip m (x m - 2atx s )\ 2 



J dx s \ip s (x s )\ 2 \ip m (x m - 2atx s )\< 



(5-11) 



of finding the signal variable x s given a measurement 
of the meter has provided x m . Moreover, the Wigner 
function of the conditional system state is now given by 



^ («Es 5 Ps I •Em ) 



1 



W{x m ) 

— OO 

x ip s (x s - |) Vm ( x m ~ 2atx s + aty) 
x £(i ra - 2atx s - aty) . (5.12) 

This Wigner function can again be expressed as the con- 
volution 

00 

P {w \x s ,p s \x m ) = J dp'W s (x s ,p s -p')Wf(x s ,p'\x m ) , 

— OO 

(5.13) 



where W s {x s ,p s ) is given by Eq. Q4.9D , and Wf{x s ,p s \x m ) 
this time reads 



W f (x s ,p s \x m ) 



1 



2wW(x m ) 

xtfj 
xip. 



dye 



x 7n 2ot ( x$ ^ 



2trf + . (5.14) 
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If we now change the variables y = 2aty in Eq. (5.14) 
and p = p'/2at in Eq. ( ^5 .13 ), we can rewrite Eq. (5.12) 
as the convolution 

oo 

P^ w \x s ,p s \x m ) = f dpW s (x s ,p s - 2atp)W f (x s ,p\x m ) 



(5.15) 

between the Wigner function of the signal state and the 
filter Wigner function 



Wf(x s ,p\x m ) 



1 



2irW(x m ) 



dye 1 



2atXe H — 
2 



xt/>* 



y 



(5.16) 



The latter can be expressed in terms of the Wigner func- 
tion of the meter [Eq.(4.14)] via the relation 



Wf(x s ,p\x Tl 



1 



W( Xr . 



-W m (x m - 2atx s ;p) . (5.17) 



Now, in contrast to Sec. VA, the filter Wigner func- 
tion ( [5.17 ) does not reduce to a delta function, and there- 
fore the Wigner function of the conditional signal state is 
not identical to the original one any more. This is indeed 
the effect of the measurement. This time, however, as we 



shall see in Sec. VI B, we can gain information about the 
signal. 



VI. METER WAVE FUNCTION 

We continue considering the meter measurement at an 
angle 9 = (p + ir/2 but discuss two extreme cases: (i) The 
meter wave function is broad compared to the signal wave 
function and (ii) the meter wave function is extremely 
narrow. In the first case we do not change the signal state 
appreciably but can only learn about the lowest moments 
of the signal distribution. In contrast, the second way of 
making a measurement destroys the state but repeated 
measurements on an ensemble of systems all prepared in 
an identical way allow us to reconstruct the signal state 
using tomographic cuts. 



A. Weak measurements 

Since ip m is broad compared to %ji a we can evaluate ip m 
at some characteristic value of x s , such as (x s ). In this 



case the conditional state, Eq. (5.9), reduces to 

oo 

|V>i c) ) = / dx s ip s (x s )\x s ) , 



and the probability 

W(x m ) = \ipm{x m + 2at(x s ))\ 2 



(6.2) 



is the original meter probability shifted by an amount 
2at(x s ). Hence, when this shift 2at(x s ) is larger than 
the width of W m (x m ) = \ip m {xm)\ •> we can learn about 
(x s ). As seen from Eq. (|Q|), in this case the state of the 
signal mode does not change appreciably. 



B. Tomographic measurements 

Optical homodyne tomography l^- [ll|j2^] is a method 
for obtaining the Wigner function (or, more gener- 
ally p9|-^l|, the matrix elements of the density opera- 
tor in some representation) of the electromagnetic field, 
preparing the field again in the same state after each 
measurement. It therefore consists of an ensemble of 
repeated measurements of one quadrature operator for 
different phases relative to the local oscillator of the ho- 
modyne detector. However, the method first employed in 
Rcf. H needs a smoothing procedure, because, in order 
to reconstruct the Wigner function one has to perform an 
integral involving the marginal probability distribution of 
homodyne measurement ||. This was indeed performed 
in Refs. |||l(j by methods which are standard in tomo- 
graphic imaging |32|. 

In the present section we show that it is possible to 
perform tomography on the meter mode to obtain infor- 
mation about the signal state. To this end, we recall 
Eq. ( pQ| ) 



W(x m ) = / dx s \ii s (x s )\ 2 \ip m (x m - 2atx s )\ 2 , (6.3) 



which gives the marginal distribution of the meter (prob- 
ability distribution of the results of the measurements of 
x m ) in the case of out of phase measurements. Let us as- 
sume that the meter wave function is extremely narrow, 
that is the meter is initially in a highly squeezed state, 
for example a squeezed vacuum |0,£), where £ = re" is 
the squeezing parameter. Then, according to Eq. (|6.3]), 
the marginal distribution W(x m ) is given by a convolu- 
tion of the modulus square of the signal wave function 
with a narrow Gaussian 



IV>' 



m y%m 1(jtx s ) | 

"l - 



l 



x exp 



Fk coshr(l 
e le tanh r 



1 



: tanh r 



— e le tanh r) 
ix^i 2(jtx s ) 



(6.4) 



Now, if the modulus r of the squeezing parameter is 
large enough, the Gaussian (6.4) approaches a delta func- 
tion in the meter and signal variables 
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\lp m (x m ~ 2atx s )\ 2 - 

and Eq. ( |6.3| ) reduces to 

oo 

1 



\2at\ 



(6.5) 



W(x m ) £* 



2cri| 
1 

2oi 



\A(x s )\ 2 s{o 



lot) 



^ V2ai/ ~ 2^ V2M) 



(6.6a) 
(6.6b) 



Hence, by measuring the probability distribution W(x m ) 
of the outcomes of the meter variable x m (for example 
via balanced homodyne detection performed on the me- 
ter field) we indirectly obtain the probability distribution 
W(x s ), up to a rescaling given by the factor 2at. How- 
ever, from Eq. ( |5.9| ) it is clear that in this case the signal 
wave function is changed, and therefore we need to pre- 
pare the signal field again in the same state after each 
measurement. This is what is usually done in quantum 
optical tomography |P-]TT|. 

The advantage of the present scheme is that we per- 
form an indirect measurement: We do not detect the 
signal mode outside the cavity (that is, we do not have 
to take the signal field outside the cavity), but we couple 
it to a meter field which is successively detected, thus 
overcoming the smearing effect introduced by the direct 
detection of the signal |||. Moreover, there is no need 
of a smoothing procedure, since we are interested in the 
marginal probability distribution W(x s ) which is directly 
related to W{x m ) through Eq. (3.6). In order to probe 
the full state of the signal field, however, we would need 
to measure the probability distribution W(x m ) for vari- 
ous values of the phase 



VII. 



NO MEASUREMENT WITHOUT A 
MEASUREMENT 



In this section we show that if a QND measurement 
performed on the signal does not alter the probability 
density of the measured observable, then the measure- 
ment process does not provide any information about 
the measured observable itself. In order to keep the pa- 
per self-contained, we prove this conclusion for the model 
considered here, but this argument holds true also in gen- 
eral, independently of the chosen model p7|. The argu- 
ment is the following: 

Let p s (0) — \ip s (0)}(ips(0)\ be the initial density ma- 
trix of the signal, and x s the measured observable, with 
x s \x s ) = x s \x s ). The initial probability density one 
would like to preserve is W°(x s ) = (x s \p s (0)\x s ), and 
we are interested in a QND measurement of x s . To this 
end, the signal is correlated to a meter which is initially 
in a certain state \ipm), and eventually a measurement 



is performed on the meter to yield the inferred measure- 
ment result x m . The measurement is then completely 
described |0| by the probability-amplitude operator 



Y {x s , x T/ 



(Xm\U(x s ,X m )\^rt 



(7.1) 



which accounts for the three stages of this measurement: 
preparation of the meter in the state \ipm), interaction 
between the meter and the signal to be measured through 
the unitary operator U(x s ,x m ) [see Eqs. (2.3) and (3.1)], 
and projection of the resulting entangled state onto the 
meter state \x m ). The QND condition [Q for a back- 
action evading measurement then reads 

[Y(x s ,x m ),x s }=0, (7.2) 

which means that x s and Y share the same eigenstates: 

Y(x s , x m )\x s ) — Y(x s , x m )\x s ) , (7.3a) 
Y f (x s ,x m )\x s ) = Y*(x s ,x m )\x s ) . (7.3b) 

After a measurement which gives the result x m , the 
system is therefore described by the density matrix 



1 



Pa 



W(x r , 



-Y(x s ,x m )p s {0)Y\x s ,x m ) 



(7.4) 



where 

W(x m ) = Tr s [F(i s ,5 m ) / 5 ;s (0)y t (x ;! ,x m )] 

dx s (x s \Y(x s ,x m )ps(0)Y^ (xs,x m )\x s ) (7.5) 



is the probability to obtain the result x m . Now, the prob- 
ability density of the measured observable after the mea- 
surement is given by 

W s {x s ) = (x s \p s \x s ) 

= w/- ^ ( x s\Y{x s ,x m )p s (0)Y f (x s ,x m )\x s ) . (7.6) 
W[x m ) 



Applying the QND condition (7.2) and ( |7.3| ) we obtain 
1 



W s (x s ) 



W(X T , 



-\Y(x s ,x m )\ 2 W^(x s ) 



(7.7) 



If we require that this probability density does not 
change due to the measurement process, W s (x s ) — 



Ws°\x), then it must be that 

\Y(x sl x m )\ 2 =W{x m ) 



(7.8) 



However, W{x m ) is not a function of x s (the eigenval- 
ues of the measured observable) and therefore also the 
eigenvalues Y(x s ,x m ) of Y(x s ,x m ) are independent of 
x s . Since the operator Y describes the measurement pro- 
cess, if its eigenvalues are independent of the eigenvalues 
of x s , the measurement obviously gives no information 
about x s , unless the measured state is an eigenstate of 
the measured observable. 
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VIII. CONCLUSIONS 

In this paper we have proposed a method to measure 
the quadrature probability distribution (or, more gen- 
erally, the full quantum state) of a single mode of the 
electromagnetic field inside a cavity. It is based on indi- 
rect homodyne measurements performed on a meter field 
which is coupled to the signal field via a QND interaction 
Hamiltonian. We have named this procedure "endoscopic 
tomography" because (i) it does not require (in contrast 
to Refs. |]-[ll]]) to take the field out of the cavity, just 
as in "quantum state endoscopy" p]], where a beam of 
two-level atoms is used as a probe; (ii) tomographic mea- 
surements performed (by balanced homodyne detection) 
on the meter mode allow us to reconstruct the marginal 
probability distribution of the signal variable or even the 
full quantum state. 

We have computed the entangled (signal-meter) state 
which arises during the evolution under the QND Hamil- 
tonian, and evaluated the conditional signal state (given 
that a measurement on the meter has provided a certain 
result). Then, we have concentrated ourselves on two 
special cases, namely, in phase and out of phase measure- 
ments. We have shown that in the first case the shape of 
the Wigner function of the signal is not changed by the 
measurement, but also that such a measurement does 
not provide any information on the signal state. In the 
second case, however, we can get information about the 
signal, but its initial state is changed due to the measure- 
ment performed on the meter: in this case, preparing the 
signal field again in the same state after each measure- 
ment, balanced homodyne detection of the meter mode 
allows the reconstruction of the original signal state. Fi- 
nally, we have given an argument according to which the 
results we have found in our model are rather general: a 
QND measurement which leaves unchanged the probabil- 
ity distribution of the system observable does not provide 
any information on the signal state. 
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APPENDIX A: DISPLACEMENT OF THE 
METER STATE 

In this appendix we calculate the state 
\i> m ) = exp[-i2atx s x m (^)][%f) m ) = D(/3)\ip m ) (Al) 



which results from the application of the operator 
D(j3) = exp[— ij3x m (4>)] on the meter state \4>m) with 
[3 = 2atx s . When we use the representation 



IV'm) = / dx m ip m (x m ;9)\x m (9)) 



(A2) 



in quadrature states \x m (6)) at the angle 9, the state 
\tp m ) reads 

oo 

\i>m) = J dx m ?p m (x m ;9)exp[-if3x m (<j))}\x m (9)) . (A3) 

— oo 

We recall the result 

cxp [-i(3x(4>)]\x m (9)) = 
"3/3 2 

—- sin [2(6» - (j))] + (3x m cos(9 - <, 



exp < — z 



x \[x m + pwa.(9-4)](6)) 
derived in Appendix M and find 



(A4) 



{V'n 



d'X rn 'Ip'm (^m 7 ^) 



3/3 2 



sin [2(9 - 4>)] + f3x m cos(6> 



x exp ^ — i 

x \[x m + p8hi(9-<l))](9)) , 



(A5) 



which after introducing the integration variable x„ 
x m + (5 sm.(9 — 4>) reads 



\lpn 



— (3 sin( 



x exp < — i 



ft 2 



sin [2(9 - <p)} + (3x m cos(0 



x \x m (9)) 



(A6) 



Hence the meter wave function gets displaced and expe- 
riences a phase shift. 



APPENDIX B: DISPLACEMENT OF A 
QUADRATURE STATE 

In this appendix we derive the relation 

e- lfim \x(9')) = exp [vp(x; /?, 9' - 9)} 

x \[x + f3$in(9' -9)](9')) (Bl) 

for the c-number [3 and the quadrature operator 



x(9) = -=(ae- w +aU ie ) 
V2 



(B2) 
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Here a and fit denote the annihilation and creation op- 
erators, respectively, with 



[a,a f ] = 1 



(B3) 



Note that according to Eq. (Bl) the action of the expo- 
nential of the quadrature operator x(9) at the angle 9 on 
a quadrature eigenstate \x{9')) at angle 9' yields, apart 
from the phase 

cp(x; 0,0-9) = — f — sin[2(0' - 9)] - 0x cos(6>' - 9) , 

(B4) 

again a quadrature eigenstate at the angle 9', but with 
the eigenvalue 



x = x + sin(0 



(B5) 



To prove Eq. (Bl) we first express the operator x{9'), 
Eq. (B2), in quadrature operators 



x{6') = -ifae^'+^e' 9 ') 
V2 



(B6) 



and 



p(9') = x(9' + tt/2) 
V2i 



ah* ') 



(B7) 



at the angle 9' . After minor algebra we find using these 
expressions the relation 

x{9) = cos(0' - 9)x(9') - sm(9' - 9)p{9') . (B8) 

The Bakcr-Hausdorff relation Ip4fl 



e A+B = e A e B e -i[A,B] 



(B9) 



for two operators A and B with [A, [A, B}} = [B, [A, B]} 
yields 

exp [-i0x(9)\ = exp[-i0cos(9' - 9)x{9')] 
x exp[i0shx(9' -0)p(9')} 



.0 2 



< exp<{ -i^-sm[2(0' - 9)} 



(BIO) 



where we have made use of \x{9'),p{9')] = i, following 
from Eqs. Q, @, and @. 

Recalling the displacement property 



\x + y) 



(Bll) 



of the momentum operator p, we find using the represen- 
tation Eq. (BIO) of the operator x(9) the expression 



cxpi ■■,.iH9)\.r(6')) : < xp<|-i^sm[2(0'-0)] 



xexpH/3cos(0'-0):z(l9')] 
x|[ir + /3sin(0'-0)](0')) , 



<-x])j- /. /.rC«)] |.r( «')} = .•x|><!-i^sin[2(e'-0)] 



exp \—%0x cos(9' 



x + 0sin(9' -9)}(9')) (B13) 



which is the result Eq. (Bl). 
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FIG. 1. Possible scheme for endoscopic tomography. Our 
goal is to measure the density operator of the signal s with- 
out taking the field out of the resonator. For this purpose, we 
couple it via a nonlinear medium of susceptibility \ to a pump 
and a meter field, p and m, respectively. The meter field is in 
a squeezed state created for example in a separate resonator 
by another nonlinear medium NL. The pump field is in a co- 
herent state of large amplitude a and phase (j>. The phase e 
of the squeezing parameter is identical to </>. Apart from the 
nonlinear coupling between signal and meter, there is also a 
linear coupling between the two. When the susceptibility a is 
equal to the product \ a i the effective interaction Hamiltonian 
for signal and meter is the product of two quadrature oper- 
ators of the two fields. In particular, the two operators are 
out of phase and their average phase is set by the pump field. 
Using part of the pump field as a local oscillator with phase 
9, we perform a balanced homodyne detection on the meter 
coupled out of the resonator. When the homodyne phase 9 
is equal to the pump phase <f> the interaction displaces the 
state of the signal field along the momentum axis and there- 
fore leaves the position distribution invariant. Moreover, in 
this case the measured meter distribution is not influenced 
by the interaction with the signal. Since we have only dis- 
placed the quantum state of the signal, we have disturbed 
it in a controlled way. Unfortunately, we have not obtained 
any information about it. When the homodyne phase 9 is 
out of phase with the pump phase <j> the interaction changes 
the signal as well as the meter field . We therefore obtain 
information about the signal field but also disturb it. Conse- 
quently, we have to reprepare all quantum states after each 
measurement. To reconstruct the quantum state of the signal 
using tomography, we record the quadrature distributions of 
the meter for all phase angles < 9 < n. In these measure- 
ments the homodyne phase 9 has to be locked to the pump 
phase 4> such that 9 = <j> + ty/2. 
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